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1 Introduction 

We recently described a new approach to the 
MEG/EEG inverse problem based on Bayesian Infer¬ 
ence [1]. Unlike other approaches to the inverse prob¬ 
lem, this approach does not result in a single “best” 
solution to the problem. Rather, it yields a proba¬ 
bility distribution of solutions upon which all sub¬ 
sequent inferences are based. Our earlier work fo¬ 
cused on the analysis of data at a single point in time. 
This work demonstrated the utility of Bayesian in¬ 
ference both for including pertinent prior information 
(anatomical location and orientation, sparseness of re¬ 
gions of activity, limitations on current strength and 
spatial correlation) and for yielding robust results in 
spite of the under-determined inverse problem. Our 
next step is to extend the spatial-only analysis to a 
spatial-temporal Bayesian inference analysis of the 
full spatial-temporal MEG/EEG data set. Preliminary 
results of this extension are presented here. 

2 Approach 

We used a parametric distributed source model for re¬ 
gions of neural activity similar to that used for the 
spatial-only analysis. Previously active regions con¬ 
sisted of voxels on cortex that were within a bounding 
sphere of some radius centered at some cortical voxel. 
Although this model was computationally efficient, it 
sometimes encompassed opposite sides of a sulcus 
into one active region. Now our active regions are 
grown from a seed voxel along tagged cortical vox¬ 
els through some number of dilations. This is a bet¬ 
ter model of our underlying picture of active regions, 
consisting of patches of contiguous cortex. Within 
each active region, neural current is now a function 
of both space and time. We also added a prior prob¬ 
ably distribution for the correlation of one time point 
with another that allows us to include the temporal 
correlation at nearby latencies. This temporal prior is 
analogous to the spatial correlation prior distribution 
we used in our earlier spatial source analysis. 

These changes have a number of effects. First, the 



Figure 1: Locations and time-courses of the three re¬ 
gions used to generated the simulated data set. The 
regions in the right hemisphere have the time-course. 


full spatial-temporal data set is now available to in¬ 
form the spatial properties of location and size of ac¬ 
tive regions. This should reduce the uncertainty of 
these parameters. Second, temporal correlation infor¬ 
mation should reduce the uncertainty of current time- 
courses within active regions as compared with not 
using temporal correlation information. Third, the 
Markov Chain Monte Carlo used to sample the many 
likely solutions now produces sample solutions that 
contain both spatial parameters of the number, loca¬ 
tion and size of active regions, and temporal time- 
courses of the total current within each active region 
as a function of time. Finally, this additional informa¬ 
tion comes at only a modest computational cost be¬ 
cause most of the computational burden is associated 
with the forward fields, which are the same for each 
latency. 

3 Results with simulated MEG data 

To test the spatial-temporal Bayesian inference ap¬ 
proach we constructed a relatively realistic simulated 
MEG data set. These simulated data were generated 
from 3 active regions shown in the three surface views 
in Figure 1. The simulated data were for a Neuromag 
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Figure 2: Waveform plot of the 122 channels of the 
simulated data set. 


122 MEG instrument. One hundred latencies were 
simulated and the time courses of the total current 
in each region are shown in the figure. The two re¬ 
gions in the right hemisphere had exactly the same 
time course. Correlated Gaussian noise was added 
with a standard deviation of 10 fT. The simulated data 
are shown in Figure 2. 

A Markov Chain Monte Carlo algorithm generated 
20,000 sample solutions, of which the final 8,000 
were used in the results presented here. Ten of these 
likely solutions are shown in Figure 3. There is less 
variability in the location and size of active regions 
than compared to the spatial-only analysis. We be¬ 
lieve this is due to the increased information content 
of the spatial-temporal data set. At least 3 regions 
were consistently present at locations and sizes that 
were consistent with the 3 true regions. Four samples 
of the total current time-course for the most posterior 
region are shown in Figure 4. Just as with the spatial- 
only analysis there is a good deal of variability, yet 
there are features of the time-courses that are consis¬ 
tent across the many samples and are therefore highly 
probable. 

A useful way to visualize the distribution of time- 
courses and to identify the probable features of the 
time-courses within each active region is to make a 
histogram of the current at each time point and to plot 
this series of histograms next to each other. This is 
shown in Figure 5 to the right for all three active re¬ 
gions. The dots indicate the true time-course values. 
Significant features can easily be identified by devi¬ 
ations of the dark band (representing the density of 
likely values) from zero. This shows the utility of 



Figure 3: Ten sample solutions from the thousands 
of solutions generated by the MCMC sampling algo¬ 
rithm. Each panel shows three views of a sample so¬ 
lution containing multiple active regions. The active 
regions are shown as maximum intensity projections 
onto the cortical surface in each view. 



Figure 4: Four sample time-courses for the most pos¬ 
terior active region. 


considering the range of likely solutions over just a 
single "best fitting" solution. In all cases the recon¬ 
structed locations, sizes and time-courses were con¬ 
sistent with the true values used to generate the simu¬ 
lated data set. 

In general, the uncertainty in the time-courses is rel¬ 
atively small. In this analysis a temporal correlation 
prior was used with a standard deviation of 3 latency 
units. We repeated the analysis but with the tempo¬ 
ral correlation removed and found this uncertainty in¬ 
creased by about a factor of 2.5. Thus, temporal cor¬ 
relation prior information may be used to reduce un¬ 
certainty in the inferred time-courses. 








Figure 5: Location and time-course distributions for 
the three active regions found to be significantly ac¬ 
tive in the analysis. The bottom row shows the lo¬ 
cations of all three regions on the left and the total 
chi-squared per degree of freedom of the data alone 
(not of the fit, so this is a measure of signal-to-noise 
of the data) on the right. 


4 Conclusion 

We have extended our previous spatial-only 
Bayesian-inference analysis [1] to a spatial-temporal 
analysis of the full spatial-temporal data set. Initial 
results with simulated data indicate that the additional 
information provided by the temporal extension sig¬ 
nificantly reduces the range or uncertainty of the 
spatial aspects of the resulting solutions. In addition, 
this additional information does not significantly 
increase the computational time of the algorithm. 
Finally, just as in the spatial-only analysis where 
we demonstrated how sampling the many spatial 
solutions could yield more robust inferences about 
the locations of active regions, the sampling of 
many temporal solutions for active regions in the 
spatial-temporal analysis also appears to yield more 
robust inferences about the time-course of regions of 
activity. 
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